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Abstract 

We develop a large-N expansion for Gutzwiller projected spin 
states. We consider valence bonds singlets, constructed by Schwinger 
bosons or fermions, which are variational ground states for quantum 
antiferromagnets. This expansion is simpler than the familiar expan- 
sions of the quantum Heisenberg model, and thus more instructive. 
The diagrammatic rules of this expansion allow us to prove certain 
identities to all orders in 1/N. We derive the on-site spin fluctuations 
sum rule for arbitrary N. We calculate the correlations of the one di- 
mensional Valence Bonds Solid states and the Gutzwiller Projected 
Fermi Gas upto order 1/N. For the bosons case, we are surprised to 
find that the mean field, the order 1/N and the exact correlations 
are simply proportional. For the fermions case, the 1/N correction 
enhances the zone edge singularity. The comparison of our leading 
order terms to known results for N=2, enhances our understanding of 
large-N approximations in general. 
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I. Introduction 



The use of large-N approximations to treat strongly interacting quantum 
systems has been very extensive in the last decade. The approach has orig- 
inated in particle physics, but has found many applications in condensed 
matter systems. Some notable examples are the SU(N) quantum Heisenberg 
model the (closely related) non linear sigma or CP^~^ models the 
Anderson and Kondo (Coqblin-Schrieffer) models H, and the two band 
Hubbard model for cuprate superconductors [0]. 

Generally speaking, the parameter N labels an internal SU(N) symmetry 
at each lattice site, (i.e. the number of "flavors" a particle can have). In most 
cases, the large-N approximation has been applied to treat spin hamiltonians, 
where the symmetry is SU(2), and N is therefore not a truly large parame- 
ter. Here lies its primary weakness, since in most cases the N > 2 models 
are not physically realizable. Nevertheless, the 1/N expansion provides an 
easy method for obtaining simple mean field theories. These have been found 
to be either surprisingly successful, or completely wrong depending on the 
system. For example: the Schwinger boson mean field theory works well 
for the quantum Heisenberg model, except for the half-odd integer antiferro- 
magnet in one dimension [|T], The latter is better described by a fermion 
large-N approximation. It is important to investigate the conditions which 
allow certain large-N generalizations to yield a "better" mean field theory for 
a particular N=2 system. 

In contrast to spin wave expansions about a broken symmetry state, 
the large-N approach can describe both ordered and disordered phases. At 
N = oo, the generating functional is dominated by its saddle point, which 
is a non interacting mean field theory with few variational parameters. The 
variational equations and the leading order correlations are in many cases 
analytically tractable. 

The corrections to the mean field theory are given by Feynman diagrams, 
where the "interactions" are mediated by RPA matrix propagators. It is 
hard in most cases to compute these diagrams even to first order in 1/N, 
which is why they have been determined only in a few select cases f^, ^. 

In this paper we shall start by deriving a new and simplified version of 
the large-N expansion suitable for evaluating spin correlations in constrained 
variational wavef unctions. These states have been used as trial ground states 
for various antiferromagnetic Heisenberg models. The exact calculation of 
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their correlations is not feasible in most cases. In one dimension, two cases 
which have been solved analytically are the Valence Bonds Solids (VBS), 
and the Gutzwiller Projected Fermi Gas (GPFG). We shall make use of 
these exact soultions in this paper. 

It is our primary purpose to study the properties of the 1/N expansion 
by using the constrained states as "toy" problems. Their 1/N expansion 
differs from that of e.g. Ref. [0] in two respects: (i) Here, the generating 
functional has no time dependence and Matsubara sums, and (ii) there is 
only one fluctuating field per site, the constraint field A, and no Hubbard- 
Stratonovich fields. Thus, we study the "pure" effects of the constraints, 
without the interactions effects of the quartic Hamiltonian. These features 
simplify the evaluation of 1/N corrections considerably. 

This paper is organized as follows. Section II introduces the valence bonds 
states. Section III introduces the Gutzwiller projected Fermi gas states. Sec- 
tion IV defines the generating functional of the spin correlation functions for 
both bosons and fermions states in a unified notation. Section V derives the 
1/N expansion of the correlation functions and describes the diagrammatic 
rules. Section VI applies the diagrammatic rules to prove three identities to 
all orders in 1/N: the absence of charge fluctuations, the sum rule for on-site 
spin fluctuations and the absence of zero momentum correlations. Section 
VII describes the results of the mean field and 1/N order spin correlations 
for the one dimensional VBS and GPFG states. The most surprising result 
is that for the VBS states of integer spin, the mean field, 0{1/N) correction, 
and the exact result for N=2 are simply proportional! Section VIII sum- 
marizes what we have learned from our approach in the context of large-N 
approximations in general. It also lists some conjectures and open questions, 
which emerge from this study. Appendices A,B, and C, fill in some technical 
details which have been used to derive certain equations in the text. 

II. Schwinger Bosons Valence Bonds States 

Schwinger bosons describe spin operators in a rotationally invariant formula- 
tion. The standard SU(2) spin operators are given by two commuting bosons 
a, h at each site as follows: 

St = ^(ajai-blbi) 
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St 

sr 



a\hi 
h\ai 



(2.1) 



The spin size s is determined by projecting the states with the Gutzwiller 
operator Vg onto the subspace which obeys the local constraints at all 
sites 



riai + riu = 2s . 



A Schwinger boson mean field wave function is defined as 



\u) = exp 



2 ^ 



Uij{alb\ — bla\) 



|0) 



(2.2) 



(2.3) 



where either determined by some mean field hamiltonian, 

or are taken as free variational parameters. It is easy to verify that due to 
the invariance of the forms ajfo] — 6|a] under global spin rotations, \u) is a 
total singlet. If we restrict to be bipartite, i.e. to connect only between 
two distinct sublattices, say A and B, then we can redefine the operators on 
sublattice B by sending 



-bj , bj 



and for i E A, j E B, Uji — 
a symmetric matrix. Under 
into 



—Uji and Uij 



u. 



J e B (2.4) 
so that u transforms into 



4f) the mean field wave function transforms 



1^) 



exp 



E 



Uij{ala] 



+ bib]) 



|0) 



(2.5) 



Eqs. ( |2.1[ - |2.5| ) can be generalized to SU(N) representations using N fiavors 
of Schwinger bosons aim, m = 0, . . . , N. In order to construct symmetric 
forms in the mean field wavefunction we again restrict ourselves to bipartite 
lattices, and define the SU(N) generalization Simm' of spin operators ( |2.1| ) as 



Si; 



i G A 
i G B 



(2.6) 



where we have generalized the SU(2) sublattice rotation 
The local constraints generalize to 



to SU(N). 



N 



m=l 



Hi 



Ns 



(2.7) 
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Ns is an integer, where s is a generalized "spin size" . 

The SU(N) generahzation of our mean field wave function ( p.5|) is 



\u) = exp 



N 



m=l 



|0). 



(2.8) 



We list some essential properties of these states in Appendix A. 

It is easy to show that for any m, m', and bond {ia,jb) we use the defini- 
tions ( p.6|) and find that 



N 



Sjmm' \ = 0. 



(2.9) 



Relations ( p.9| ) show that \u) is globally SU(N) invariant, and is therefore a 
singlet of total spin. 

We shall restrict ourselves to translationally invariant states, which in 
Fourier representation are given by 



|-u) = exp 



N 



keBZ m=l 



|0), 



(2.10) 



where BZ is the first Brillouin zone, Wk = J2j e^'^-'woj, aim = A/" ^ J2j ^^^''(^]mi 
and N is the number of lattice sites. We also define Skmrn' = Z]j Sjmm'- 



II. 1 The Gutzwiller Projection 



The mean field state ( |2.8| , |2.1CI| ) includes different spin sizes at each site. In 
order to construct a bona-fide state of spins s, we must project out all other 
spin sizes using the Gutzwiller projector 



\ll)s =Vs\u) 



(2.11) 



which enforces the constraints ( |2.7|) . By expanding the exponential in (|2.8| ) 
and applying the Gutzwiller projection, we obtain 



1/1 ^ 
1^). = 2^«..E 



|0) 



m=l 



5 



, , 1 X! n ''-'^i ( X! ^im^^m ] 



|0) 

(2.12) 



where i^b = \ J^Ns is the total number of bonds in the projected state and Ca 
labels the different configurations ofuf, bonds on the lattice, where exactly Ns 
bonds eminate from every site. In Fig. 1 we depict several configurations for 
various {uij} and values of Ns. \u)s which is a sum over such configurations, 
is called a "Valence Bonds State" . 

Since all Simm' commute with the constraint, \'u)s is also rotationally 
invariant and a total singlet. If all the bond parameters are non negative, 
Uij > 0, the wavefunction (|2.12|) satisfies the Marshall sign criterion [§. We 
recall that the ground state of any bipartite Heisenberg antiferromagnet must 
be a total singlet and obey Marshall's theorem. 

The correlation function in the Gutzwiller projected state is defined as 

^ {h j) ~ {Simm' Sjm'm) (2.13) 

and its Fourier transform is 

5--'(k) = ^5'""^'(0,j)exp(zkj) 

j 

= {S -kmm' S-km' m) : (2-14) 

where for any operator O we denote 

{O) ={u\rsOVs\u)/{u\Vs\u). (2.15) 

II. 1.1 Resonating Valence Bonds (s = |) 

Several special cases have received particular attention in the literature. For 
N = 2, s = 1/2 and nearest neighbors \u)s is a superposition of all 

dimer configurations. In one dimension there are only two such (Majumdar- 
Ghosh) configurations []TU|, which have an exponentially small overlap for 



large lattices. The spin correlation in one dimer state vanishes beyond the 
nearest neighbor range. In two and higher dimensions, s = | states are sums 
over many valence bonds configurations, which were denoted as "Resonating 
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Valence Bonds" (RVB) by Anderson [jrT|. One configuration in the square 



lattice RVB state is depicted in Fig. l.a . The RVB state was proposed 
by Anderson []TT| and others as trial ground states for frustrated quantum 



antiferromagnets and high superconductors. The number of dimer config- 
urations grows exponentially with A/", and the overlap between different con- 
figurations is finite. The computation of the spin correlations in the dimer 
and longer range RVB state was carried out numerically by Liang, Doucot 



and Anderson |T2| using bipartite bonds of various range. They found that 
the RVB states with Uij ~ I/t^" have long range order for a < 5. 

II. 1.2 Valence Bonds Solids, Integer s 

Affleck, Kennedy, Lieb and Tasaki [|13| (AKLT) found a class of extended 
SU(2) Heisenberg models for which the exact ground states are Valence Bond 
Solid (VBS) states, given for SU(N) by 

/ N 

vbs\ TT / \ ^ t „t 



i^''^^ = n E 




im jm 

{ij) \m=l 

M = Ns/z, (2.16) 

where (ij) denotes nearest neighbor bonds and z is the lattice coordination 
number. The condition that M must be an integer restricts the size of the 
spin and lattice for which such states can be defined (see Figs, l.b, l.c). 
For example, the SU(2) model in one dimension allows s = 1,2,3.... On 
the square lattice, only even spins s = 2, 4, 6 . . . are allowed. The correlation 
function S+-{i,j) = for one dimension has 

been calculated for all s by Arovas, Auerbach and Haldane (AAH) to be 



S^-{^,J) = i-iy-''^-^^^ exp[-K,\J-^\] 
= log(l + 2/s), 

s'-w - lis+i) ^'-;;^f, . (2.17) 

3 1 + cos{k) + 

The real space correlations decay as a pure exponential, with 1/ks as the 
correlation length. is a "spin liquid" ground state of the kind that 



was predicted by Haldane |T^ using the large-s, Non Linear Sigma model 
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analysis of the Heisenberg antiferromagnet. A AH also found a Haldane gap 
ISj] in its single mode excitation spectrum. The correlations of VBS states 
on higher dimensional lattices are those of a classical logarithmic Heisenberg 
model, at temperature T = z/ s. This implies that for large enough s, the VB 
states in three dimensions have long range Neel order. The calculation of Eq. 
(|2.17| ) was performed in the SU(2) coherent states basis. The generalization 



of this calculation to > 2 has not yet been achieved. In the following, we 
shall apply the large-N expansion to this problem. 

It may be verified that the one dimensional Schwinger boson state \u''^')s 
with 

uf =5(,,) , uf' =2cos{k) , (2.18) 
is dominated by the VBS state in the limit of infinite lattice size M: 

(M!)5-^" \u''^')s = \'^'''")s + 1^'), (2.19) 

where (\1/'|\I'') ~ for some c > 1. The exponentially small corrections 
are of non uniform valence bonds configurations, where some bonds have 
higher powers of a^a^ than others. Consequently, we expect that in the 
thermodynamic limit N' —>■ oo the spin correlation function in the state 
|'U^*'^)s is given also by ( |2.17| ). 



III. Gutzwiller Projected Fermi Gas 

In this section we introduce another important family of variational states 
using fermions rather than Schwinger bosons, to represent the SU(N) spin 
operators 

The local constraint on the fermion occupation is 

N 

Him = rii = Ns, (3.2) 

m=l 

where A^s is an integer, which by the Pauli principle must be less than or 
equal to A^. Using the fermion operators, one can construct a global SU(N) 
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singlet by the following state 




(3.3) 



where kp is the Fermi momentum which is chosen to include Ns states per 
site in the Fermi volume. Vs Gutzwiller projects onto the subspace which 

Uk are variational parameters, 
in an exponential form as follows: 

1 



satisfies the constraint (|3 
It is possible to write 



\u)s = Vs exp 



Y^u-kOikp - |k|) 51 4m 



t 



-km 



|0), 



(3.4) 



where Wk = — W-k- In real space ( |3.4| ) are analogous to the Schwinger boson 



states, defined in (|2.11| ), ( p.8|) , where Uij = Af ^ J2\]i\<kF '^k^ 



jk(j-j) 



U 



—Uij. A well known case is the Gutzwiller Projected Fermi Gas for s 



[i.e. a half-filled Brillouin zone 

■f 



=Pi exp 



^sign(k)^(A;^-|k|) 



a} a} 

"km "-km 



|0). (3.5) 



In real space |v|/9p/9^ contains long range bonds Wjj, \i — j| >> 1. Since 
the bonds are not bipartite, it does not satisfy the Marshall sign criterion. 
This state is deduced from the mean field theory of Baskaran, Zou and An- 
derson 1^ for the Heisenberg antiferromagnet. In one dimension, |vI>9p/9^ 
for SU(2) was found to be the exact ground state of the Haldane-Shastry 
hamiltonian |T^, whose interactions fall off as the second inverse power of 
distance. This state has correlations, similar to that of the ground state of 
the nearest neighbor Heisenberg model. Haldane has also shown that the 
Haldane-Shastry hamiltonian and the nearest neighbor Heisenberg model 
share similar gapless excitation spectra |[T7|] . 

Gebhard and Vollhardt |[T^ have calculated the correlation function of 
one dimensional I^sp-^^') for N=2 



log 1- 



which in real space decay asymptotically as a power law 

ZTT J ~ ^ 

where Si(x) is the sine integral function. 



(3.6) 



(3.7) 
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IV. Correlations and the Generating Func- 
tional 



The spin correlations of\u)s can be derived from a generating functional. The 
generating functionals for the Valence Bonds states ( p. 11 ) and the Gutzwiller 



Projected Fermi Gas states (|3.4|) are formally very similar, and given by 



Z[j] = { exp 



imm' ^im^im' 



(4.1) 



where are either bosons or fermions, jimm' are the source currents, rji = 1 
for fermions and for bosons 



1 ieA 
-1 ieB 



(4.2) 



which takes care of the sublattice rotation of the SU(N) spins (|2.6| ). The 
functional derivatives of Z determine the spin correlation functions. It is 
sufficient to use symmetric source matrices jimm' = jim'm- Hence, jimm' 
and jim'm o,re not to he considered as independent, but should be varied 
simultaneously when differentiating Z[j]. 

The following relations can be directly verified from ([4.1|) 



5 log Z 



Sji 



ij^im} 



(4.3) 



j=0 



which is a direct consequence of the SU(N) symmetry in ( |2.8| ) and the con- 
straint ( p.7|) . The two-point spin correlation function Eq. ( p.l3| ) is given 

by 

1 „ . 6^Z 



2 ^mm' ^jimm'^jjmm 

Additional terms generated by the differentiation in 

{Sim:^m' Sjm^m') 0. 



(4.4) 

j=Q 

must vanish, since 
(4.5) 



Eq. (|4.5| ) follows from the rotational invariance of the wave function. It is 
easy to verify, that Sim^m'\u)s and Sjm'^m\u)s are eigenstates of the operator 
J2i{Simm — Sim'm') with eigenvalues +2 and —2 respectively. Therefore these 



two states are orthogonal and (|4.5|) follows 
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For m 7^ m', the correlation function S"^^"^' {i, j) is an SU(N) general- 
ization of the usual SU(2) spin correlation function S~^~{i,j) = (S^S^). In 
rotationally invariant states this function is related to the correlations of the 
other SU(2) spin components by 



=2{S^S]) =-(S.-S,) 



(4.6) 



The evaluations of Z and S"^"^'{i,j) of Eqs. ( [4.1| , [4.4|) are complicated 
because of the Gutzwiller projector. If it were absent, we could easily cal- 
culate Z as a matrix element of an exponential bilinear operator as done in 
Appendix A. In order to proceed, we must choose a convenient representation 
for the projector. The projector can be represented as a limit of a strongly 
interacting density matrix. 



Vs = limexp 



1 



(4.7) 



Keeping e finite will help to control infrared divergences in subsequent dia- 
grammatic calculations. The matrix elements of ( |4.7|) are hard to evaluate 
in its present form. Using an auxiliary constraint variable Aj at every site, 



we transform ( [4.1|) to an integral 

v. = lun£mexp[E(- 



X^i+i\i{ni - Ns) 



PA ^ n 




(4.8) 



Now we can write the generating functional as 



lim / VX {u\ exp a^{iX+j)a 



1^) 



X e 



where we denote the matrices 



(4.9) 



A Xi5 a' 5 jYiYn' 1 



3 Vijimm' ^ii' 



(4.10) 
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In ( [4.9| ) we have used the commutation of the spins with the density operator, 







(4.11) 



to combine the exponentials of the source terms and the projector. Now we 
use Appendix A, to evaluate (|4.9|) as 



Z[j] = lim|PA exp{NS[\,j]) 
C 



2N 



Tr,^log(l-CMV^+^Me*^+^) 



IS 



(4.12) 



where C = +1 (~1) ^oi bosons (fermions) and is the hermitian conjugate 
of the matrix u. 

The correlation function is given by Eq. 



^"'"'(1,2) 



1 .i„„z-' fvx (n- ^'^'^-^l 



2 — 6mm' 



where S[X] = S[\,j] 



j=0 



+ N' 



6j\mm' 6j2mm' 

,6S[X,j] 6S[X,j] 



6jlmm' 6j2mm' 



exp {NS[X]) 



j=0 



(4.13) 



V. The Diagrammatic Expansion of Z{j) 

The multidimensional integration over VX is equivalent to the difficult com- 
binatorical problem of evaluating the correlations in the valence bonds state, 
e.g. ( |2.12| ). This is clearly seen by expanding the action in powers of e*'^* 
and integrating with the weights e~'^^^^\ The integrals reduce to products 
of (5-functions, which select the terms with Ns powers of Uij for any given 
site i. Summation over these terms is very cumbersome in general. 

In ([4.12| , |4rT3| ) the parameter was scaled out of the action S. We shall 
evaluate the A integrals by a saddle point expansion which is controlled by 
the largeness of A^. The functional S[X] is expanded as a Taylor series about 
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its minimum A; the coefficients of expansion are independent of N. Since 
the wave function is translationally invariant we shall search in the space 
of uniform saddle points Aj = A. A is found by requiring that the linear 
variations vanish. We define e^'^ = u. The saddle point or "mean field" 
equation is 

1SS[X] 



i SXj 



u ■u}u{l — (u ffu) 



33 



0, 



(5.1) 



which determines u[u,S]. We see that ( |5.1| ) implies that u is real, i.e. the 
integration paths of the variables Aj have to be analytically continued to 
cross the imaginary axis at A = —ihi{u). The term [ ]jj in ( p.l| ) is the 
average number of particles of every flavor in the (unconstrained) state \uu) 
(see Appendix A). Thus, the mean field equation yields that the constraint is 
satisfied on the average. If we define the unprojected state using uu instead 
of u, the mean field equation (|5.1| ) would be satisfied with A = 0. Hence we 
shall use that convention for u and set A = 0. 
We now expand S, in powers of A to obtain 

s[\] =5(°) -ly-c(2) 

where iS™* includes only third and higher order terms: 

oo -n 

Cint _ Q{n) \ \ 

^ - 1^ -T'^n,...,i„ ^ii ■ ■ ■ -^in- 



ij AjAj + iS™*, 



n=3 



(5.2) 



(5.3) 



g{n) depend parametrically on u. Diagrammatically they are depicted 
as thick circles with n A vertices denoted by wiggly lines (see Fig. 2 ). Later 
we shall obtain explicit expressions for iS^"-* in form of loo'ps constructed of 
Greens functions. 

We also expand the pre-exponential functions in ([4.13 ) 



1 



Sji mm' 



'Vi 



oo -n 



3=0 



N 



(n+l) 



Smm'Vl ^ i 



■ ■ Ai + S 



2 ^mm' ^ilmm'^32mm' 



3=0 



N 

V1V2 
N 



(n+l) 



n=0 



n! 



A, 



^ q{n+2) . . 

2^ -T'^l,2,ji,...,i„ ^ii ■ ■■^i 



n=0 



(5.4) 
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where sums over repeated indices are assumed and we denoted S[^^ = s. Here 
we use (|4.12|) to relate the derivatives with respect to A to those with respect 
to j. Diagrammatically, we denote the current vertices (at points 1,2) by 
dashed hues as shown in Fig. 2. 

The correlations functions can be evaluated by inserting the expansions 
(lOD and (U) into (|413|) 



where 
^^(1,2) 



// oo -n 



(5.5) 



n=0 
oo -n 



X 



(n) 



L=0 



\n=?t 
oo -n 



N 



oo -n 



X 



(n) 



L=0 



L\ 



■■■Xi 

Ail ■ • ■ Ai., 
il ■ ■ ■ Aj 

A,;, ■ ■ ■ Ai 



exp 



\n=0 "'■ / \n=0 "'■ 



\n=3 



exp ( -y4?AfcA/^ 



We disregard the contribution of the constant NS^^^ in the action. 

The integrals in (|5.6|) are sums of multidimensional gaussian integrals. 

The gaussian integrations contract all A fields in pairs A^A; bringing down a 
propagator for each pair given by 

c« =^a; =-i(s'^))-;, (5.7) 

where the minus sign is due to the factor of i, which accompanies every A field. 
The propagator ( |5.7|) is depicted as a wavyline connecting two A vertices in 
Fig. 2. One has to sum over all A contractions. The disconnected parts of 
the diagrams serve to cancel the factor of Z~'^, leaving us with the diagrams 
which are connected to one or the other current vertices (A linked cluster 
theorem). Thus calculating any particular diagram involves multiplying loops 
and propagators D, and summing over internal lattice points. There are 
internal loops, created by powers of 5™*, and external loops, coming from an 
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expansion of the preexponential functions ( |5.4| ). Internal loops must have at 
least three A vertices; external loops have current vertices and might also have 
arbitrary (including zero) number of A vertices. The order of any particular 
diagram is given by 



where L is the number of internal loops (or the power L of the sums in ( p.6| 
and P is the number of propagators (half the number of A fields in ( p.6| 
After grouping all the diagrams at each order in 1/N we obtain the series 



(5.8) 



oo 



^"^^"'(1,2) = ^Ar-P^'"'"'(p)(l,2). (5.9) 

Similar rules govern the calculation of higher correlation functions. One 
has to sum over all possible ways of distributing the current vertices on the 
external loops. Within each loop the m, m'-indices of the external currents 
must be equal to the indices of other external currents, to allow non zero 
values of the trace. 



VI. Identities To All Orders 

The diagrammatic expansion of the 1/N series has special structure which 
allow us to obtain exact identities to all orders in 1/N. A key feature is that 
the propagator D of the constraint field A is non other than the inverse of 
the square part of the action S^'^\ As a result we shall show that the local 
constraints are exactly enforced to each order in the expansion, i.e. there 
are no contributions to charge fluctuations when all terms of the same order 
are considered. In addition we shall prove a sum rule for the on site spin 
fluctuations for arbitrary and the absence of zero momentum off-diagonal 
spin correlations. 

VI. 1 Absence of Charge Fluctuations 

Here we shall demonstrate that the constraint is imposed at each order in 
1/A^. In other words, due to the Gutzwiller projection the density fluctu- 
ations vanish identically after all diagrams of a given order are summed, 
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yielding 

{niA) = Ns (A) (6.1) 

for an arbitrary operator A. It is instructive to see how (|6.1|) is derived by 
the diagrammatic expansion. A current of ni belongs to some external loop 
^(n+i)^ n > 0. Let us first consider all the contributions with n > 1. We 
define a "tail" of a diagram as the combination of a propagator attached 
in series to a loop S*^^-* which has the operator rii on its other vertex. All 
diagrams can be separated into two classes: ones with a tail, and ones without 
a tail. It is easy to identify for each diagram without a tail say R{ni,A), a 
counterterm ^(^i,^) by attaching a tail to the rii vertex (See Fig. 3). By 
(|5.7D the two are of the same order p (they have the same number of loops 



minus propagators) and they cancel precisely 

R{n,,A) = NY^SflDuMnuA) 

kl 

= -R{ni,A). (6.2) 

Thus, at any order p, the counterterms cancel the connected charge fluctua- 
tion diagrams one by one. The only terms in the expansion of {riiA) which 
survive are the disconnected contributions with rii on the loop 5^^). An im- 
portant property of the diagram rules is the absence of the counterterm to 
the S^^^ loop. Such counterterm would involve an internal loop S^^\ which is 
not allowed by the rules of our expansion. Thus (|6.1|) follows from iS'-^^ = s. 
Q.E.D. 



VI. 2 Sum Rule for the On-Site Spin Fluctuations 

First we shall show that for SU(N) invariant states there is a relation between 
two types of correlation functions: diagonal 5'™'" and off-diagonal S"^^"^ . 
Both functions do not depend on the particular values of m, m' since the wave 
functions are SU(N) symmetric. Due to the Smm' coefficient in ( 
is equal to . Let S^"^ be the connected part of the diagonal correlation 
function. Any diagram which contributes to the off-diagonal correlations also 
appears in S^"^, i.e. 

if i?" G S""^"'' , then G S^"". (6.3) 

However, in S^"^ its contribution is partially cancelled by the counterterms, 
which are given by tails insertion. It may be seen, that the sum of any 



5.5), S"^^' 
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diagram G S'^^"^' and its counterterms, obtained by all possible ways of 
insertion of tails, is equal to (1 — l/iV)i?". 

In addition to diagrams of the type i?" and their counterterms, yet an- 
other contribution to S"™™ is given by the diagrams R^, which have the two 
current vertices on different loops S^"^^^^^ and S^'^^~^^\ ni^2 > 2, and by their 
counterterms . But the diagrams R^ must exactly cancel with their coun- 
terterms! This is verified by adding a tail to one of the loops and seeing 
that the counterterm is of the same order in 1/N due to the additional m- 
summation for the loop which became an internal loop. This proves that 
all the diagrams which do not cancel are of the type i?" and the following 
important identity holds: 

= (i-^) ^'"^"'(^-j)- (6.4) 

Now we can calculate the onsite spin fluctuations. For SU(2) we are 
familiar with the "spin square" operator S^, which when projected to the s 
sector yields a diagonal matrix of elements s{s + 1). For larger A^, its natural 
generalization is 

= 'Pns ^ ^ ^imm'^im'm 
mm' 

"Pns ^ ^ ij^im'iX ~^ C^im) C^mm'^im') 
mm' 

= {n^s{1 + Cs)-CN) Vns, (6.5) 

where in the fermion case we made use of equality = (C = — !)• On 
the other hand, using ( |6.4|) we can write the onsite fluctuations as 

N{N -I) + N {1-^)^8'^^'^' +Ns^. (6.6) 
Using ( |6.5|) we find the desired identity 

5-^-'(^,0 =]^^(1 + C^)- (6-7) 
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For SU(2) this formula reduces to the known values: 



/q+q-\ _ / + 1) Schwinger bosons , . 

{^^^^) - I i fermions (s=l/2) ^^'^^ 

In translationally invariant cases ( |6.7| ) gives a sum rule, which is useful for 
checking the diagrammatic calculations at each order of 1/N seperately. In 
momentum space, the diagrams of order must obey 

=(l.)\_QPs(i + (:s), p = 0,l,2,... (6.9) 



k 



VI. 3 Absence of zero momentum correlations 

The last identity is a consequence of the singlet nature of the wavefunction 
which implies that 5'k=om^m'|^)s = 0. By using Eq. ( p.l4| ) we obtain from 



that 

5™^"^'(k = 0) = , (6.10) 
which holds of course at each and every order in the 1/A^ series. 



VII. Calculations of leading orders 

In this section we calculate the spin correlation functions 5""^™' using the 
1/A^ expansion. We start with an explicit evaluation of the loops 5^"^ and 
propagator D. In the cases of interest u is hermitian. It is useful to define 
the following Greens functions 

= C^/'^uu (^±1 - C^/\uy\ (7.1) 

We also introduce the matrix A = 5ij{e^'^^ — 1) and express the action ( [4.12| ) 

as 

S[\] = 5(°)-|5:Trlog(l-t.,A)-^.5:\-^5:A,^ 

7=± i i 

5W = -^5:Trlog(l-7C'/'^n). (7.2) 

^ 7=± 
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is a constant which we shall disregard. 
Expanding the logarithm in (7^) and using to cancel the linear term 
we obtain 

A oo 2 

'^[A] = |E -TT{u,Ar--j:^l (7.3) 

where X)' denotes that terms linear in A are excluded. By equating terms of 

the same order in A in Eqs. ( |5.2| ) and ( |7.3| ) we can relate the loops S^"'\ n > 2 

to traces over Greens functions u^. Diagrammatically, we denote the Greens 

functions by thin sold lines. A closed loop of Greens functions denotes a trace 

over lattice and 7 indices. For a n A vertices iS^"-*, there are contributions 

from diagrams with 1 < m < n Greens functions, since the function A yields 

all powers of A fields at the same point. Due to Eq. (|7.4|) below, the loop 
5(1) 

may be denoted diagrammatically in the same manner - as a closed loop 
of one Greens function with one vertex. 

The translational invariance of u makes it easier to work in the momen- 
tum representation. The linear action, or the mean field equation (|5.1| ), is 
explicitly given by 

I + =jf^ = ^'-'^ 

The quadratic action in (|5.2|) is given by 

2 c 

^ 7 

= 9E«7^J%•^ + '^^J^ + ^^Je^ 

7 

5^'^ (q) = X E E «7k «7k+q + 5 + 

^■/V k 7=± 

k .til (7 - c^/%^k) (7 - c^/^^«k+q) + ^ + ^ ' ^^-"^ 

and the propagator is 

Z}(q) = - (Ar5(2)(q))^\ (7.6) 



19 



An important property of this expansion is that the lowest order (mean field) 
correlation function 5*"^^"^ (°)(q) is simply related to the quadratic part of 
the action: 

om^m' (0) / N 

J5(2)(q + 7f), Schwin erer bosons «n 
~ i5(2)(q), fermions ^ ' ^ 

where vf = (vr, tt, . . .) for a cubic lattice. The diagrams for 5''"^™'' (q) are 
shown in Fig. 4. At this point we note that -D(q) is singular for e ^ since 
^m^m'(o)(^q = 0) vanishes by ( |610D . This causes diagrams which involve one 
or more propagators to diverge as 1/e. A check on the correctness of the 
calculation is that these "infrared" divergences must exactly cancel between 
different diagrams to yield a finite result for lim^^o 'S'(^''(l, 2) for each order 
p separately. We shall come back to this point in our summary. 

The 1/N corrections for 5"*"^™ (k) are given by the diagrams of Fig. 5. 
Solid lines represent factors of u^{k). Each vertex conserves momentum, and 
indices 7 = +, — of the solid lines. We must sum over internal momenta and 
7. Diagrams with external currents 1,2 at the same point denote an overall 
factor of 612. Due to the cancellation mechanism described in section VI. 1, 
the third and fourth diagrams in the bottom row of Fig. 5, cancel against the 
fifth diagram. Thus by the same mechanism, there is complete cancellation 
between the last four diagrams in the second row. We shall describe the 
calculations of the remaining diagrams for the Scwhinger boson and fermion 
cases separately, and defer technical details to Appendices B and C. 

VII. 1 Valence Bonds Solid Correlations 

The mean field equation Eq. ( [7. 41 ) for the Valence Bonds Solid state 
is 

dk 4u2cos2(fc) 



2n 1- 4m2cos2(/c) 



whose solution is 



(7.8) 



Js{s + 2) 
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By (|73) and (|7J) we obtain 

^m^m'{o)^j^^ ^ _2^2 y- r (k cos(g) cos(fc + q) ^ ^ 

^' ^t^^J-TT 27r (7 - 2Mcos(g))(7 + 2ucos(fc + g)) 

(7.10) 

The integral is performed by introducing a new variable z = e*^ which trans- 
forms the integration over q into an integration along a unit circle |z| = 1 in 
the complex 2;-plane. Using { \l.9\) yields 



Sm^-'^ik) = {s + l) ' If ^^ , . (7.11) 

1 + cos(A;) + 

This result is very surprising, since it is just proportional to the exact result 
for = 2 as found by Ref. and given in Eq. ( p. 171) . The factor 2/3 



between (|2.17|) and the mean field result ( [7.11| ) is consistent with the factor 
N/ {N + 1) between the mean field on site fluctuations and the exact sum rule, 
Eq. ( |6.7|) . This suggests that perhaps the exact dependence of 5""^™ [k] 
is given by this simple multiplicative factor. Fortunately, we are able to 
calculate the 0(1/A^) corrections and check this proposition at least to the 
next leading order in 1/A^. This calculation is described by the diagrams 
of Fig. 5, which involve exchanges of one or two propagators D. The sum 
of all diagrams was evaluated analytically using the symbolic manipulation 
program Mathematica [^. The result, derived in Appendix B, is 

^m^m'{l)(^^ = -^5™^™'W(fc). (7.12) 

This result confirms the above hypothesis, but is far from obvious! In fact, 
the separate 1/N diagrams have infra-red divergences of order due to the 
diverging propagators at momentum tt. In addition, the separate diagrams 
have different correlation lengths than the mean field function, but these 
effects somehow cancel by summing all the terms of order 1/N, leaving us 



with Eq. ( 7.12 ). It is highly tempting to conjecture that the same relation 
holds to all orders, i.e. that 

S'^^'^'^P\k) = (^-iy^-^'"'W(A;) (?) (7.13) 

which will sum up to the simple relation 

S^^^^'ik) = s'^^'^'^^^k) (?) (7.14) 
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For N=2, we have already seen that this conjecture is correct. But, as we 
shall discuss in section VIII, the underlying reason for this relation is still a 
mystery. 



VII. 2 Fermions GPFG Correlations 



The mean field equation ( [7. 41) for the one dimensional GPFG state |v|/sp/9^ 



5D is given by the integral 

■u^ r dk 



——^ / -e{k,-\k\) -s = 0, (7.15) 
with A;^ = 7r/2 and s = 1/2. Eq. ( f^lSD yields 



n{s) = (7.16) 



TVS 



which implies that for kp — * 7r/2 and s = l/2,'u— >oo. To enable us to 
calculate the correlations for the GPFG, we keep u finite by holding a Fermi 
level kp a. little above it/2. We shall take the limit kp —>■ 7r/2 only at the end 
of our calculations. This divergence simplifies the calculations considerably. 



because Eq. ( |7.1|) yields a simple limit 



TT 



hm^u^{q) = -e^--\q\j. (7.17) 
The lowest order correlation function is evaluated using ( [7. 7] ) , ( |7.5| ) and ( [7.17] ) 



fcF^f ' ' 2 J-^27r V2 "7 V2 

\k\ 



9^ ^^-^S) 

ZTT 



When comparing ( |7.18| ) to Gebhard and VoUhardt's result, Eq. (P^), we 



find that the two expressions agree very well in the small k limit, where they 
vanish with the same linear coefficient, but they deviate at larger k as shown 
in Fig. 6. Eq. ( |3.6| ) diverges logarithmically near k = n, while the mean field 
result (|7.18|) has merely a discontinuity in its derivative. This translates to a 
difference in the asymptotic power law decay in real space, between l/\i — j\ 
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of Eq. (^) and - jf of (^AE^). There is a factor of A^/(A^ - 1) = 2 
between their sum rules as required by Eq. 

In Appendix C, we calculate the diagrams for the GPFG state. We 
obtain the result: 



(7.19) 



In Fig. 6 we compare the functions S""^""' ^^^k) , S'""^""' ^^\k) + 
and the exact result Eq. ( |3.6| ) for N=2. We see that the correction 
improves the mean field approximation considerably near the zone boundary, 
where its derivative diverges logarithmically. In real space we obtain for 
separations r = \j — i\, 



1 


'\k\ 










(, - SI) „ 




N 









[7 + log(7rr) — Ci(7rr)] 



where 7 = 0.577... is the Euler constant, and Ci(x) = — J!^ dt cos(t)/t 
vanishes for large x. In ( [7.20|) we find that the 1/N correction enhances the 



long distance correlations from r ^ to r ^ log(r). 



VIII. Summary and Discussion 

In this paper we have introduced a large-N expansion for the correlation func- 
tions of Gutzwiller projected states. We have discovered several properties 
and sum rules which hold to all orders in By explicitly calculating the 

mean field and 1/N corrections for particular Schwinger boson and fermion 
states, we can check the validity of this approach against exact results for 
N=2. We shall conclude by discussing what we believe we have learned from 
our results, and what still needs to be illuminated by further investigations 
and insight. 

VIII. 1 What We Understand 

1. The effects of the Gutzwiller projector can be expanded systematically 
in terms of 1/A^ diagrams. Each diagram with L loops and P propaga- 
tors is of order (1/A^)^^~^). The loops and propagators are determined 
by the mean field (saddle pont) equation. 
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2. 



The local charge fluctuations are suppressed at each order by countert- 
erms, which have the tail structure depicted in Fig. 3. 



3. For any A^, the onsite spin fluctuations are given by 




N + C 



s{l + Cs) 



(8.1) 



where C = 1 bosons (fermions). 

4. Each diagram can diverge due to the divergence of the propagators 
(see discussion after Eq. ( [7.7|) . However, the sum of all diagrams of the 
1/N order is finite. We conclude that in general, large-N expansions are 
prone to such intermediate divergences, due to the "hardness" of the 
constraints (or lack of "self-interaction" for the A fields). The lesson to 
be learned is that results which are based on any subset of diagrams, 
or on partial resummation schemes, are highly suspect. 

VIII.2 What We Do Not Understand 

1. For the VBS states, the mean field, 0{1/N) and the exact N = 2 
correlations are simply proportional. We conjecture that the higher 
order terms behave in the same manner, i.e. 



For Schwinger bosons, we know that this relation holds for the on-site 
sum rule (|8.1| ), but its validity for all /c is a surprise. We can recall 
however that similar surprises have been found in other large-N calcu- 
lations, both with bosons and with fermions, where mean field results 
differ from the exact result by the factor N/ {N + (). For example: The 
mean field susceptibilities of the s = | ferromagnet in one dimension, 
and antiferromagnet in two dimensions [|I| (both for N = 2) are off 
by a factor of 2/3 = N{N + 1). Also, the Wilson ratios of the Kondo 
impurity model and the s = ^ Heisenberg antiferromagnetic chain 
m are 2 = N/ [N — 1). It would therefore be very useful to understand 
this relation in order to correct the mean field approximation for other 
problems. The apparent simplicity of this correction factor may have 




N 



S^^-^'i^){k) (?) 



(8.2) 



A^+ 1 
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its origin in some group theoretical relation between the saddle point 
approximation and exact integrals over Haar measures ||21||. 



The above discussion indicates that for these systems the 1/N expan- 
sion is not just an asymptotic series but a convergent, well behaved 
expansion. On the other hand we are faced with the apparent failure 
of the boson large-N theory for the Valence Bonds Solid at s = | in one 
dimension. The series yields exponentially decaying correlations, 
while the correct state (the nearest neighbor dimers state) has van- 
ishing correlations beyond nearest neighbor separations. We therefore 
strongly suspect that there is an essential singularity in the expansion 
of the form 

^ [1 + exp{i27cNs)] (8.3) 

which cannot be obtained at any order in the expansion. Such a factor 
distinguishes between integer and half odd integer spins for N = 2. 
This is similar to the famous topological Berry's phase, or "^-term" , of 
the continuom theory of half odd integer Heisenberg antiferromagnets 
in one dimension. This term must be added to the Schwinger boson 
mean field Lagrangian to obtain the correct ground state degeneracies 



3. We note that the Fermion large-N approximation is quite successful 
for the s = i GPFG state in one dimension. The corrections 
enhance the long distance correlations from to r~^log(r). It would 
be interesting to find out how the full 1/N series modifies the power 
law to for N = 2. 

We recall that the fermion mean field theory for the spin half Heisenberg 
chain |TB[ is successful in reproducing the Fermi-liquid features of the 



exact solution . Here we have found another empirical evidence that 
the fermionic approach is better than the bosonic approach for s = \ 
antiferromagnets in one dimension. In two dimensions, the relative 
advantage of the fermionic versus bosonic large-N approach is not clear. 
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Appendix A: The Properties of Mean Field 

States 



The major goal of this appendix is the derivation of Eq. ( ^4.12| ). We will 
define here the mean field states \u) by 



\u) = exp 



2 ^ --ij^i^j 



ij 



|0), 



(Al) 



where Oj are either bosons or fermions, satisfying the usual commutate rela- 
tions 

(A2) 



a,- a 



(ajtti = ajaj — (a^al = 



( = +1 (—1) for bosons (fermions) and matrix u satisfies the symmetry 



condition = ( 



u [u 



Uji). We have dropped the m indices for simplicity. 



The states |-u) have the following important properties: 
1. The overlap of mean field states is given by 

-- 

{u\v) = det(l — (u^v) ^ . 



(A3) 



2. Let A is an arbitrary operator, the product of any number of creation 
and annihilation operators. Then the extended Wick's theorem holds: 
the normalized matrix element of A, {u\A\v) / {u\v) , is equal to the sum 
of all possible completly contracted products of creation and anihilation 
operators (with the usual sign for fermions), in which each contraction 
involves two operators. A contraction of di, d2 is the normalized ex- 
pectation value {u\did2\v) / {u\v) . 



3. The contractions are given by 



{u\alaj\v) 
{u\v} 



U 



{v\aW.\u) 
{v\u} 

{u\v) 



(A4) 
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The proof of these properties may be found in ||2^ (see also ||2 



Following we will write the whole set of creation and annihilation 
operators as a 2A/'- dimensional vector (A/" is the lattice size): 



7 = {a, a^} = {ai, ...,a^f,a\, 
with commutation relations 

TiTj Cljli Pij) 

where p is the 2M x 2M matrix 



We have = so p 



-1 



-(p. Following Balian and Brezin 









J = exp 


^7^7 


= exp 



(A5) 
(A6) 

(A7) 
we define 

(A8) 



where the 27V x 27V matrix R satisfies the symmetry condition = (R (in 
this appendix we will call such matrices symmetric) and 27V x 27V matrices 



T = 



Tu Ti2 



T21 T22 



(A9) 



It is shown in p2[, that the matrices ([A9|) faithfully represent the second 
quantized operators ([ASf ), i.e. 



J{Ri)J{R2) = J{R) - 
We will now prove the inverse statement, 

• J{R^)J{R2) = J{R) , 



^pRi^pR2 _ gP-R 



(AlO) 
(All) 



that is to say, the representation of operators JT" by matrices T is isomorphic. 
To prove the (|A11|) , we use the Baker - Campbell - Hausdorff formula, which 
gives 
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e=^e=^ = exp I ^ " " " [Q., ■■■]]), (A12) 

, n=l li,...,l„=l,2 
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where we denoted ci,2 = |7-Ri,27 and are some constants, which we 

do not need to know exphcitly. On the other hand, since e^^^e''^^ = e''^, the 
same formula gives 



R = P^Y. Oil,^...,ln[pRln[pRln-l---[pRl2,pRh]- 
n = l li,...,ln = l,2 

Correspondingly, ( [All| ) will be proven, if one can show that 
, 1 



(A13) 



Let us denote 

An = p'^[pRi„[pRi„_, ■ ■ ■ [pRi^,pRi, 



^p-'[pRi„[pRi„_, ■ ■ ■ [pRi,,pRi,] ■ ■ ■]]!. (AM) 



(A15) 



so that the r.h.s. of (|A14| ) is equal to ^^An^. We will first prove by induction, 
that An is symmetric, = (An- For = 1 it is correct; then, for An+i we 
have 

An+l = P'^ [pRu+i , pAn] = Rln+ipAn - AnpRl„+, , (Al6) 

and since -Rz„+i is symmetric and p is antisymmetric, p^ = —(p, the symme- 
try of An+i follows from the symmetry of An- 

Now we will prove by induction the relation ( [A14| ). For = 1 it is trivially 
correct; then, if it is correct for n, we have for n + 1: 



(A17) 



Using symmetry of -R/„^i and A^ and the commutation relation for 7 (|A6|) it 
is straightforward to show, that r.h.s. of ( [A17| ) is equal to 

ifp^' [pRi„+^ , pMi = \ip~^ [pRin+i [pRin ■ ■ ■ [pRh , pRh] ■ ■ (AI8) 

which completes the proof of ( [Allj ). Note, that we have also shown, that the 
product of two operators of the type ( [A^ ) is another operator of the same 
type, represented by the matrix R of (|A13|) . 

Now we shall use this "multiplication rule" to calculate ( |4.9|) . The expec- 
tation value is of the form 



M 



{u \ exp 



\u) = (OlJiexp 



J2\0), 



(A19) 
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where J\ = exp 



and J2 = exp 

8]) and 17 is A/" x A/" matrix. We can transform exp 
by writing 

1./ c^^ \ . 



are operators of the type 
also into the form 



exp 



exp 



rri ^ 

-— ir ?7 
2 ' 



exp 



2'^ \ V 



7 



exp 



C rri 

-- ir 77 
2 ' 



(A20) 



where jTa is of the type ([A8D, so we have for M 

(o|JiJ3:^2|o) 



M 



exp 



det 



-— ir r? 
2 ' 



(0|JiJ3J^2|0). 



(A21) 



The T matrices, corresponding to J\, J2 and jTs are 

' 



Ti 



To 



exp 



exp 



exp 



1 

-C 

1 

-C 

1 

-C 





' 

u ^ 

Cr)^ 

1) 













:) 



1 u \ 
1 j ' 



e-'' 



(A22) 



Thus, T, which corresponds to J" = JiJ^J2 is 



T = T1T3T2 



e ' M 



(A23) 



It is shown in ||2^ , that for every operator J ([A8| ) and corresponding matrix 
T the following relation holds: 

(0|J|0) = [detT22]-i (A24) 

Using Eqs. ( [A24|) , ( |A23| ) and ([A21|) we will obtain for the case of symmetric 
(in usual sense) 17, fj^ = fj, the final relation 



{u \ exp 



\u) 



det (l-Cuh^ 

c 



ue 



exp 



-^Trln(l -Cw^e^ ue^ 



(A25) 
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which yields the expression ( [4.12| ) for the generating functional ( [4.9| ). 



Appendix B: Calculation of ( |7.12| ) 

Here we explicitly calculate the order diagrams for the Schwinger boson 
case. The diagrams are depicted in Fig. 5. We use the integration variable 
z = e^^ instead of the momentum k. For example, for the nearest neighbor 
bonds problem if^'^ of ( p.l8| ), we have 



u''^'{z) =z+-. (Bl) 
z 

The conservation of momentum at every vertex is equivalent to the rule that 
the product of all z's entering a vertex is equal to unity. Each sum over k 
[ p5|] is replaced in the thermodynamic limit by a contour integration over z 
on the unit circle 

lim U-^yPk = — r dkFk ^ —i —Fiz) 

N^oo ^ 2nJ-7r 27liJ\z\=l z ^ ' 

= Y.^es[F{zi)lz,] (B2) 

i 

For the quadratic part of the action (|7. 5| ) the integral is 

S^^{z) = ^ E ^ / '^-^u,{z')u,{zz') + s + e\ (B3) 

where 

u±{z) = uu{z) [±1 - uu{z)r^ (B4) 

and for the valence bond case u{z) is given by (pl|) and u by ( [7.9| ). 

Since trigonometric integrands F^ are replaced by rational functions F{z), 
it is easy to determine their poles Zi (including a pole at z = 0), and their 
residues at these poles. The sum over residues in ( p3D yields 

S^'\z) = is+ l)n + e^ (B5) 

[Z — K)[l — ZK) 

where 
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The propagator D{z) is equal to — (^NS^'^^z)"^ . The diagrams are gener- 
ated by the rules of Section 3, and shown in Fig. 5. As an example, the 
integrations of the vertex diagram (see Fig. 7) are 

1 1 f dz' f dz" 

2 5^ J^~^ J~jJ7' "^^^'^ u^{zz') u^{zz'z") u^{z'z") D{z") . (B7) 

Each diagram in Fig. 5. diverges as 1/e. However, the divergences 
cancel in the sum, and the overall 1/N result is finite for e — > 0. The sim- 
plicity of the residues method allowed us to use the symbolic manipulation 



program "Mathematica" |^ to perform the integrations analytically on the 
computer. The program identifies the poles and residues of the rational func- 
tions. Intermediate expressions, especially for the diagram Fig. 7, involved 
upto hundreds of terms. Expanding these terms and finding common de- 
nominators became too cumbersome for manual calculations, and therefore 
automating this process was essential. The result of this calculation is given 
by Eq.dni). 



Appendix C: Calculation of ( |7.19| ) 

Here we derive the order 1/N correlations of the Gutzwiller Projected Fermi 
Gas state, Eq.( [7.19D . The quadratic part of the action is given by ( [7. 181 ), so 



by (|7.6| ), we have a diverging propagator at /c = 0. We control this divergence 
by letting kp > 7r/2. We denote 

^-l^f>0. (CI) 

The role of 6 is similar to that of e - it regulates the divergence of the 
propagator. Thus we obtain 

5{2)(fc) = ^±M, (C2) 



so that 



27r 



We first note that the last two diagrams of the top row in Fig. 5, cancel since 
they differ by one line, which yields a factor of —1. For the same reason, 
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the sum of the first two diagrams is equal to twice the contribution of the 
second diagram. The remaining contributions to S"^^"^ ^^\k) are depicted 
(including combinatorial factors) in Fig. 8. By reflection symmery, we can 
restrict ourselves to A; > 0. 

The first diagram of Fig. 8 is given by 



Ri{k) 



D{q)P{k + q) 



dq 
2^' 



(C4) 



where P{k) is the polarization bubble. By Fig. 4, for k in the first Brillouin 
zone 

P{k) = S^^\k)-sJ^-\ 

which yields for 



(C5) 



Ri{k) 



N 



k\ 7T — k k k 2k 
1 log — log - + 1 

n I n n TV 



(C6) 



The contribution of the second diagram of Fig. 8 is 



2R,ik) = 2 f ^D{q) f ^e{^-\p\)e{^-\p + q\)e{^-\p + q + k\). 

(C7) 

We denote the integral over p in (|C7| ) as (27r) ^A{q), where 



which yields 

2R,{k) = i- [ 

JM \_ TV J-TT —q 

The last diagram of Fig. 8 is 

dq 



TT + g if— 7r<g< —k 

n-k if-A;<g<0 
{-K - q-k)d{'K - q-k) if < g < vr 



(C8) 



1 f-^ TT + q 



1 -K-k 



I X '^^ , 

TT J-k — q Tf Jo 



Rsik) = - r p-D{q) r ^eC--\p\)e{^-\p + k\) 



1 [■^-'^ 71 - q - k 
5 + q 



dq 



(C9) 



xO{l-\p + q\)e{^-\p + k + q\). (CIO) 
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The integral over p is equal to {2n) ^{tc — k — \q\) 6 {ji — k — |g|), so ( |C10| ) 
cancels the last term of (|C9|) , while first two terms of ( PoD yield 



2R^{k) + R^{k) = - 



' 5 k ^ 5 ^ k' 
log log - + 1 

-n Tx k IT . 



(Cll) 



By adding (^11]) and we obtain Eq. (^A^). 
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Figure Captions: 



Fig. 1: Graphical representation of valence bonds configurations, which 
contribute to (a) Resonating Valence Bonds States (RVB) on the square 
lattice, (b) Valence Bonds Solid on the square lattice, (c) Valence Bonds 
Solid on the chain. 

Fig. 2: Diagrammatic representation of the 1/N expansion for the 
correlation function S{1,2). Thick circles denote loops, wiggly lines are 
A-propagators, and dashed lines are external currents. 

Fig. 3: Diagrammatic representation of cancellation of charge fluctuations 
by counterterms, Eq. (|6.2|). 

Fig. 4: The diagrams contributing to the loop S^'^\ u± are mean field 
Greens functions, defined in Eq. (|7.1| ). 

Fig. 5: The 1/N corrections to S'"^^'"'(k). See discussion after Eq. (|7.7|) . 
Fig. 6: The N=2 spin correlations of the Guztwiller Projected Fermi Gas 



in one dimension. The exact result (solid line, from Ref. [0) diverges at 
k = 71. The mean field (MP, short dashes) result has discontinuos derivatives 
at /c = 0, TT, while the sum up to order 1/N (long dashed line) has a diverging 
derivative at = vr (see Eq. (|7.19|) ). 

Fig. 7: Assignments of integration variables in Eq. (|B7]). 

Fig. 8: Diagrams contributing to the correlations of the GPPG (Ap- 
pendix C). 
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